The role of interactions in spin-polarised atomic Fermi gases at unitarity 
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We study the zero temperature properties of a trapped polarized Fermi gas at unitarity by assum- 
ing phase separation between an unpolarized superfluid and a polarized normal phase. The effects 
of the interaction are accounted for using the formalism of quasi-particles to build up the equation of 
state of the normal phase with the Monte Carlo results for the relevant parameters. Our predictions 
for the Chandrasekhar-Clogston limit of critical polarization and for the density profiles, including 
the density jump at the interface, are confirmed with excellent accuracy by the recent experimental 
results at MIT. The role of interaction on the radial width of the minority component, on the gap 
of spectral functions and on the spin oscillations in the normal phase is also discussed. Our analysis 
points out the Fermi liquid nature of these strongly interacting spin polarized configurations. 
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I. INTRODUCTION 

The behaviour of spin-polarized Fermi systems has 
been a fascinating subject of research for over fifty years 
since it can lead to the interplay of superfluidity and 
magnetism. In particular a Fermi superfluid is resistant 
to spin polarisation and so one can naturally ask what 
happens when we attempt to polarize it. 

One such example occurs in a superconducting metal 
when we apply a magnetic field. Under certain condi- 
tions the coupling to the orbital motion (responsible for 
the Meissner effect) is negligible and the important effect 
is the coupling to the electron spins. The field can lower 
the energy of the spin-polarized normal state and, if it is 
strong enough, make the normal state energetically more 
favourable than the superconducting spin-singlet state. 
The value of the field at which this transition takes place 
is known as the Chandrasekhar-Clogston limit [l| and, 
in a BCS superconductor, it requires that the field (or, 
in neutral systems, the chemical potential difference be- 
tween spin up and down) be larger than A/y/2 where A is 
the gap. Crucially, this estimate assumes that the change 
in energy of the normal state due to polarization is only 
kinetic in origin and neglects changes in the interaction 
energy. 

However, if the system is strongly interacting, the value 
of the Chandrasekhar-Clogston field can also depend on 
the interactions in the normal state and we must accu- 
rately take them into account if we wish to study nor- 
mal/superfluid coexistence. 

Exactly this situation has arisen in recent experiments 
that have been carried out in the strongly interact- 
ing unitary limit of two-component atomic Fermi gases 
@) S 0, S @1- Such gases can be polarized leading to 
a state most naturally described as phase separated be- 
tween a normal and a superfluid component. The ques- 



tion then revolves around the energetics of such a system. 
Most of the previous theoretical work has concentrated 
on the nature of possible superfluid phases and assumed 
that interactions in the normal state are unimportant 
0- 

In the present work we analyse in detail the role of 
interactions in establishing the phase diagram of the uni- 
tary spin-polarised system as well as in other relevant 
phenomena as the density profiles of trapped gases and 
the frequency shifts of the collective oscillations. We 
will develop the Fermi liquid theory of the normal state 
which was introduced in Q and show that not only its 
predictions for the critical polarization (Chandrasekhar- 
Clogston limit), but also the detailed structure of the 
density profiles (including the density jump at the in- 
terface between the superfluid and normal phases) agree 
very well with the experimental data obtained at MIT 
[H, The same theory provides explicit predictions 
for the frequencies of the spin excitations in the normal 
phase. The experimental measurement of these frequen- 
cies would provide a direct measurement of the interac- 
tion coupling constants characterizing the normal phase 
itself. 

We begin in Section [II] by reviewing the theory of the 
normal state and we study the resulting nature and prop- 
erties of the superfluid/normal phase transition. Then in 
Section IHII we review the key points of the MIT experi- 
mental situation concentrating on those which are crucial 
in understanding the phase diagram at very low tempera- 
ture. Next we make a detailed comparison of theory and 
experiment, clarifying the features that interactions are 
responsible for. In Section HVl we study the spin modes 
in the normal phase and make predictions for their fre- 
quency shifts due to interactions in the unitarity limit. 
Finally in Section [V] we draw our conclusions. 
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II. INTERACTIONS IN THE SUPERFLUID 
AND NORMAL PHASE 

We concentrate on T = configurations in the unitar- 
ity limit, when the scattering length, between the spm-f 
and spin- J. species, is infinite and the system, in many 
respects, behaves like a strongly interacting fluid. Impor- 
tant features accessible experimentally are related to the 
polarization of the Fermi gas which permits to explore 
different quantum phases. In particular recent experi- 
ments have revealed the occurrence of phase separation 
between a unpolarized superfluid and a polarized normal 
gas 0]. 

In what follows we shall consider only two phases: the 
unpolarized superfluid and the partially polarized normal 
phase. If not otherwise specified, we will call ng = ngf = 
rigx the density in the superfluid phase and simply n-f 
and nj, the normal state spin-| and spin- J, densities. 

The equation of state for a homogeneous unpolarized 
superfluid at unitarity is simply given by 



N s 



&^(6^ns) 2 / 3 ^6 S (n s ), (1) 
5 2m 



where Ns is the number of atoms in the superfluid phase 
and the universal parameter £g = 0.42 has been cal- 
culated employing Quantum Monte Carlo simulations 

A convenient way to build up the equation of state for 
the normal state is to consider a dilute mixture of spin-j. 
atoms immersed in a non-interacting gas of spin-t atoms. 
When the concentration x = ni/n-\ is small the energy 
of the system can be written in the form 



— ^ = -cft 1-AaH x b ' 6 

iV T 5 T V m* 

- 3 ( \ 



Bx' 



(2) 



where JVf is the number of spin-| atoms and e^t = 
S 2 /2m(67r 2 n-f ) 2 / 3 is the Fermi energy of the spin-t gas. 
The first term in Eq. ([2|) corresponds to the energy 
per particle of the non-interacting gas. The linear term 
in x gives the single-particle energy of the spin-down 
particles, while the a; 5 / 3 term represents the quantum 
pressure of a Fermi gas of quasi-particles with an effec- 
tive mass to* . Eventually, the last term includes the 
effect of interactions between quasi-particles. The val- 
ues of the coefficients entering in Eq. @ have been 
the object of various studies using different many-body 
appoaches[l, [ll|, [TH, [ID, [H, EH- In the present work we 
use the most recent values A = 0.99(1), m* = 1.09(2) and 
B = 0.14 calculated in [l5[, using Fixed- Node Monte- 
Carlo techniques. It is worth remarking that, while Eq. 
@ is thought as an expansion of the normal state energy 
for small concentration, it agrees very well with Monte- 
Carlo calculations also for large values of x Q . 

We assume that the Fermi gas is confined by an 
isotropic harmonic potential V(r) — muj 2 r 2 /2 and that 



a non-zero polarization 



P 



(3) 



can give rise to a phase separation between a unpolarized 
superfluid core and an external partially polarized normal 
shell. We define Rs the surface which separates the two 
phases and the external surface of the (fully polarized) 
normal part. In the local density approximation the free 
energy reads 

E = 2 [ dr( es (n s (r)) - M ° + V(r))n s (r) 

Jr<R s 



dr(3/5e,Fje(x(r))nj (r) 

R s <r<R^ 

+V(r)( ni (r) + n T (r)) - M °n T (r) - ^(r)) (4) 

where /ij are the chemical potentials of the spin-|, 
spin- 1 component respectively and ^ = (/x° + fi®)/2 is 
the superfluid chemical potential. 

The equilibrium is found by minimizing the energy 
with respect to the densities of the superfluid and the 
normal part as well as to the border Rs between the two 
phases. By varying the densities we find the LDA ex- 
pressions for the superfluid and for the normal densities 



^ = £s^(67T 2 n s ) 2/3 + nr), 



(5) 
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From Eq.([6]) we find that the external radius of the nor- 
mal part is simply given by R 2 = 2/j® /(muj 2 ). 

We are dealing with a polarized system and the same 
atoms are constituents of both the superfluid and the 
normal phase, thus the number of particles in the two 
different phase is not fixed a priori but is determined by 
the chemical and the mechanical equilibrium. The lat- 
ter is obtained by varying the energy functional Eq. (j4| 
with respect to Rs-, yielding the equal pressure condition 
between the two phases reads 



" 2 de S ' 
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r=R s 



From Eqs. ([5][5]), we obtain an implicit equation for the 
critical concentration at the border, namely 



e(x(R s )) + -(l-x(R s ))e'(x(R s )) 

- (2^ 5 (e(x(R s ))) 2/5 = 0. 



(9) 
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Since in the superfluid phase the densities of the spin- 
| and spin- J. component are equal, the occurrence of a 
solution of Eq. (|9|) with x(Rs) ^ 1 reveals the exis- 
tence of density jumps at the interface. Using the above- 
mentioned values for A, m* and B we find x(Rs) = 0.44, 
which also corresponds to the critical concentration for 
a homogeneous system to start nucleating the superfluid 
phase[8|. We find that the jump is significant only in the 
minority component being 



ns(Rs) 
ni(R s ) 
ns(Rs) 



2£s 



e(x(R s )) 



3/5 



x(R s 



2gs 
e(x(R s )) 



1.02, 

3/5 



0.45. 



(10) 
(11) 



In addition to the critical concentration Eqs. ©-([7]) 
can be used to calculate the density profiles which lead us 
to the important prediction that when the polarization 
is larger than the critical value P c = 77% the superfluid 
core disappears and only the normal phase is present in 
the trap [||. 

A deeper insight on the role of interaction can be ob- 
tained in the high polarization case, i.e., iVf ^> iVj,. Start- 
ing from Eq. ([6]) and ([7]), to the leading order in x, we 
get 



2m' 



Mt = ^(6A) 2/3 + nr), 



lM = ^(6-V) 2/3 +nr)(i 
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(12) 



(13) 



yielding the non-interacting value (j2 — {(SN^Y^huj for 
the chemical potential of the majority component. The 
above equations reveal that both the majority and the 
minority components have an ideal Fermi gas profile, the 
latter being described by a renormalized mass m* and 
feeling a renormalized external potential. The radius of 
the minority component is quenched by the interaction 
to the value 



R l - R °l 



^ ( , 3 A 

— ll + -A 

m \ 5 



-1/4 



(14) 



where R\ = (ASNi) 1 ^ ^h/(mu) is the Thomas-Fermi ra- 
dius of the ideal Fermi gas. These results can be also un- 
derstood by introducing the effective single quasi-particle 
Hamitonian 



■fM+|^+nr)(i + -i 



(15) 



for the minority component. 

The above formalism allows to characterize, in terms of 
the interaction parameter A, the energy threshold of the 
spectral function J2 n \ < n\c\ c^ iP |0 > \ 2 S(uj + E n — Eq) 
for the pseudospin flip of the minority component to a 
third (hyperfine) state |3 >. This function is relevant 



for photoemission spectroscopy [16[ and RF experiments 
fl7l ]. where it enters integrated over the momentum p 
(for recent theoretical calculations for highly polarized 
samples see [l8l [l9|). 

In the absence of final state interactions the threshold 
in the center of the trap is given by 



A(threshold) = \a^ + 2- (l - 

5 1 2m V m* 



(16) 



which gives in principle a method to measure directly the 
effective mass. 

Other measurable quantities which can give informa- 
tion on the parameters are the frequencies of the spin 
collective oscillations of the normal phase. Once again 
in the highly polarized case they can be easily calculated 
using the Hamiltonian Eq. (jT5]) as we discuss in Sec. IIVI 



III. COMPARISON WITH EXPERIMENT 

As the theory developed in the previous section is 
based on a local density approximation it is important 
to focus on experiments where such an approach, which 
ignores details of the interface at the microscopic scale, 
is applicable. This favours the comparison with the MIT 
experiments of [H, H| carried out with large values of par- 
ticle number N. Due to strong anisotropy and smaller 
values of N, the Rice experiments reported in [5j, [6( are 
instead expected to be much more sensitive to surface 
effects [20j and they will not be considered here. 

The main results carried out at MIT at the lowest tem- 
peratures at unitarity can be summarized as follows: 

1) There is a critical polarization P c ~ 75% above 
which there is no evidence of superfluidity. This value has 
been established by studying vortices 0] and using the 
ramp method into the BEC limit to measure the number 
of condensed molecules. Furthermore the disappearance 
of the condensate coincides with the disappearance of a 
central region where the spin-| and spin-| densities are 
equal Q. 

2) At values of the total polarization P < P c three 
regions can be distinguished in the trap: a central su- 
perfluid core in which both densities are equal, a first 
shell where both spin-| and spin-| atoms are present but 
with different densities, an exterior shell where only spin- 
| atoms are present. 

3) For P > P c only two regions are left: an interior 
partially polarized core and the exterior fully polarized 
shell. 

4) At the interface between the superfluid and the par- 
tially polarized phases there is a large jump of around 
40% in the density of the down atoms whereas the dis- 
continuity in the up atoms, if any, is too small to be 
detected. In going from the partially polarized to the 
fully polarized phase the densities are continuous. In the 
spirit of LDA we therefore identify the former transition 
as first order due to the abrupt density change, whereas 
the latter one is identified as second order. 
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All the four above-mentioned facts agree very well with 
the theory for the normal phase exposed in the previous 
section. 

We proceed by comparing the theoretical density pro- 
files with the available data[3|,|4( . In Fig. ((!]) we compare 




FIG. 1: Column density profiles for the majority, the minority 
component and for their difference: MIT data vs. our theory 
(curves obtained by the method described in the text) for 
P = 44%. 



our theory with the experimental data for the column 
density for P — 44% and in Fig. @ we show the com- 
parison with the double integrated density difference 

/oo 
dpp(ni(z,p) - ni(z,p)), (17) 
-oo 

for a range of polarisations both above and below P c ~ 
77%. As it can be observed, the matching is excellent in 
all cases. 

The agreement between theory and experiment is also 
revealed in Fig. ([3]) where we analyse the reconstructed 
three dimensional density: the jump in the nj, density 
is evident (see also 21]). We emphasize that, since the 
data are normalized to the non-interacting gas, the only 
input needed in Figs. (HJ and ([3]) is the polarization P. In 
the case of Fig. (|5J), where the experimental data are not 
given in dimcnsionlcss units, we have used the measured 
central (z — 0) density as a fitting parameter. Below the 
critical polarization the doubly integrated density pro- 
files exhibits a typical plateu whose existence is a direct 
consequence of the LDA in the presence of an unpolarized 
superfluid core [22J. 

We stress that the good agreement is ensured by the 
proper inclusion of interactions in the normal state. This 
is not just a quantitative question since assuming that 
the normal state is noninteracting can lead to unphysi- 
cally high values of polarisation for the Chandrasekhar- 
Clogston limit. This is the case of the extensively em- 



doubly integrated density 
[10 a atoms/cm] 




FIG. 2: Doubly integrated density profile difference (rid), data 
(black dots) and our theory (red lines) for different polarisa- 
tions, above and below P c , starting with the lowest curve and 
movino- nre P = 58%. 73%. 8(1% 92% 




FIG. 3: Density profiles for a polarization P = 44%. Theory: 
solid black line (dashed red line) is the spin-| (spin- J.) density. 
Experiment: the black (red) line is the spin- j (spin- J.) density 
as reported in Q. The density jump in the J. component is 
clearly visible. 



ployed Bogoliubov - de Gennes theory. Although the en- 
ergy of the superfluid even at unitarity is quite close to 
the Monte Carlo one, for the normal phase such a theory 
includes the interactions only through the pairing terms 
in the superfluid phase yielding a value close to 100% for 
the critical polarization [7j, a value definitely ruled out by 
the MIT experiments. 



IV. SPIN MODES 

The excellent agreement with experiments discussed 
in the previous section strongly supports the basic pic- 
ture underlying this paper, i.e. the assumption that, at 
unitarity, a phase separation takes place between a super- 
fluid unpolarized core and an external polarized normal 
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phase. It is useful to further exploit the relevant features 
of this normal phase. In this section we will discuss the 
dynamic behavior and we will show, in particular, that 
the study of the out of phase oscillations of the spin-up 
and spin-down components can provide a further crucial 
test of the role of interactions in the normal phase. 

A first crucial point to discuss is the role of collisions. 
Collisions are quenched by Pauli blocking at very low 
temperature, so that at sufficiently small temperature 
the small amplitude oscillations in the normal phase are 
described by collisionless dynamics [23|. In the oppo- 
site collisional regime of relatively higher temperatures 
the dynamics is instead governed by the equations of hy- 
drodynamics and the out of phase oscillations are over- 
damped since the spin current is not conserved in this 
case. Notice that the possibility of studying the collision- 
less behaviour in the normal phase of a unitary Fermi gas 
at very low temperature is a unique opportunity provided 
by spin polarized Fermi gases. In the non polarized case 
the gas is in fact always superfluid at such temperatures. 

We will consider the normal phase above the 
Chandrasekhar-Clogston limit and in particular we will 
focus on the simplest case of a highly polarized gas where 
the number of particles in the spin-down component is 
much smaller than in the spin- up component (iVj, <C Nf). 
In this limit the collective oscillations can be classified 
into two cathegories: the '"in phase"' oscillations, where 
the motion is basically dominated by the majority com- 
ponent, and an '"out of phase"' oscillations where the 
minority component moves in the trap, the majority one 
being practically at rest. The most important example 
of the first class of oscillations is the dipole center of 
mass motion. This motion is not affected by the in- 
teraction because of the translational invariance of the 
force and occurs exactly at the frequency ljd = where 
i = x,y, z. A second important "'in phase"' oscillation is 
the quadrupole mode. This mode is affected by interac- 
tion effects except in the deep N± <§; limit where its 
frequency approaches the value ljq = 2ujj_ (for simplic- 
ity we consider only the radial mode in an axi-symmetric 
trap). Notice that this frequency differs from the hydro- 
dynamic value \/2lo±_ [24j]. 

Let us now discuss the "out of phase" (spin) oscilla- 
tions. According to the discussion of Sec. [Til in the 
Ni -C Af limit these oscillations are described by the 
single quasi-particle Hamiltonian Eq. lTSl) . 

The frequency of the spin oscillations are hence easily 
calculated. The spin dipole frequency is given by @] 

while the spin radial quadrupole frequency is given by 

The above results point out that interactions affect the 
spin frequencies in a sizable way by increasing the corre- 



sponding values in the absence of interaction by a factor 
~ 1.23. Their measurement, together with the determi- 
nation of the radii and/or the spectral function thresh- 
old, would consequently provide unique information on 
the separate value of the relevant interaction parameters 
to* and A. In this respect it is worth noticing that the 
quantities to*/to and 1 + (3/5)A enter in Eq. (|14[) and 
in Eqs. (|18H19p via a different combination. Notice also 
that differently from the radius (fT4f and the threshold 
(frnj), the spin frequencies (|18[) and (fH?)) are independent 
on the number of atoms. 

The above results can be generalized to include larger 
values of N±/N^, through the proper inclusion of the in- 
teraction term between the minority and majority species 
in the energy functional introduced in Sec. |TTJ 

An important question is how these spin modes can be 
excited experimentally. It is worth noticing that an exter- 
nal perturbation coupled to the total density n = n-\ + n± 
will never excite the spin dipole mode, but only the center 
of mass motion which is a pure density oscillation, insen- 
sitive to two-body forces (Kohn's theorem). Since at the 
high magnetic fields characterizing the resonant regime 
the magnetic coupling is practically the same for the two 
hyperfine states of 6 Li used in the experiments, only 
an optical coupling with a laser field suitable detuned 
with respect to the internal atomic frequencies of the 
two hyperfine states can consequently excite the dipole 
spin mode. The situation is different for the quadrupole 
oscillation. In fact in this case, due to the absence of 
the analog of the Kohn's theorem, the spin excitation is 
not completely decoupled from the density probe and a 
quadrupole perturbation of the form Xiyi will result 
in the excitation not only of the "in phase" mode, but 
also of the spin mode (|19[) with a relative weight propor- 
tional to < r\ >i Ni/ < r\ > T A T . 



V. CONCLUSIONS 

By assuming phase separation between a superfluid 
and a polarized normal gas, we have provided an accurate 
picture of the zero temperature behavior of a trapped po- 
larized Fermi gas at unitarity. Using the recent results 
of ab initio Monte Carlo calculations for the relevant 
interaction parameters, in both the superfluid and the 
normal phase, we have found excellent agreement with 
experiments both concerning the value of critical polar- 
ization (Chandrasekhar-Clogston limit) and the shape of 
the density profiles, including the density discontinuity 
at the interface between the two phases. The agreement 
confirms that the basic structure of these novel config- 
urations is now well understood. In particular it shows 
that the proper inclusion of the interactions is a crucial 
requirement for a quantitatively and qualitatively correct 
description of the polarized phase which emerges like a 
new example of strongly interacting normal Fermi liquid, 
consisiting of a weakly interacting gas of quasi-particles 
and remaining normal even at the lowest temperatures. 
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The excitation of the spin oscillations in the collision- 
less regime (see discussion in Sec. HVj) . the measurement 
of the radial width of the minority component (see Eq . 
(fT4")l ). the availability of photoemission spectroscopy [16[ 
and RF measurements |l7( and the behaviour of the gas 
under the constraint of an adiabatic rotation [25| are ex- 
pected to provide further insight on the role of the inter- 
actions in these polarized Fermi gases. 

We finally point out that there are a number of is- 
sues related to the border between the superfluid and 
the normal phase. In particular when mesoscopic effects 
become important (small number of atoms, strong trap 



deformation) the proper description of the interface, be- 
yond LDA, can play an important role in characterizing 
the phase separation between the superfluid and the nor- 
mal component. 
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